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A fundamental connection between superconductivity and quantum spin fluctuations in underdoped 
cuprates, is revealed. A variational calculation shows that Cooper pair hopping strongly reduces the 
local magnetization mo. This effect pertains to recent neutron scattering and muon spin rotation 
measurements in which mo varies weakly with hole doping in the poorly conducting regime, but 
drops precipitously above the onset of superconductivity. 
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When holes are introduced into the copper oxide planes 
of high T c cuprates, spin and charge correlations change 
dramatically. The local magnetization too, measured 
by ^iSR and elastic neutron scattering ^ on e.g. 
La2_ x Sr a; Cu04, reveals a qualitative difference between 
the insulating and superconducting phases: mo is rather 
insensitive to doping in the poorly conducting regime 
< x < 0.06, but drops precipitously above the on- 
set of superconductivity at x > 0.06, becoming unde- 
tectable at optimal doping x « 0.15. Theoretically, holes 
can cause dilution and frustration in the Heisenberg an- 
tiferromagnet, which create spin textures: either random 
("spin glass") or with ordering wavevector away from 
(7r,7r) (sometimes called "stripes") ||. However, the ap- 
parent reduction of local magnetization by the onset of 
superonductivity, is a novel and poorly understood ef- 
fect. Theory must go beyond purely magnetic models, 
and involve the superconducting degrees of freedom. 

We find that this problem is amenable to a variational 
approach, using hole-doped Resonating Valence Bonds 
(RVB) states, originally suggested by Anderson for the 
spin correlations of high T c cuprates ■ 

These RVB states are excellent trial states for doped 
Mott insulators, with large Hubbard repulsion U : 

(i) Configurations with doubly occupied sites are ex- 
cluded. 

(ii) Marshall's sign criterion for the magnetic energy ||] 
is satisfied, and Heisenberg antiferromagnetism at zero 
doping is accurately recovered. 

(iii) For doped systems, spin and charge correlations are 
parameterized independently, without explicit spin nor 
gauge symmetry breaking. 

These are important advantages over commonly used 
Spin Density Wave, Hartree-Fock and BCS wavefunc- 
tions for the antiferromagnet, metal and superconducting 
phases respectively. RVB states permit an unbiased de- 
termination of ground state spin and charge correlations 
appropriate for the cuprates. 

A phenomenological low energy effective Hamiltonian 
is used, with two major components: Heisenberg inter- 
action for spins, and single or Cooper pair hoppingng 
kinetic energy for fcrmion holes. 

Our key results are as follows: 
(i) For the magnetic energy alone, the local magneti- 



zation too is weakly dependent on doping concentration. 
This holds independently of inter-hole correlations for ei- 
ther randomly localized or extended states, 
(ii) In contrast to (i), too is strongly reduced by the ki- 
netic energy of Cooper pair hopping, which correlates the 
reduction of too with the rise of superconducting stiffness, 
and hence |^| the transition temperature T c . 

Our results agree with the experimentally reported cor- 
relation between too and T c |QB]. This relation appears 
to be weakly dependent on the precise hole density. A 
brief discussion concludes the paper. 

The Wavef unctions: The hole-doped RVB states are 
compactly defined by 



^[u, v; x] = Vg(x)iP[u, v 



exp 




ip[u,v] 



V G {x) 



Where a\ , b\ and // are Schwinger bosons and hole 
fermions respectively |10|, i = 1, . . .L is a site index on 
the square lattice, and Vcix) is the Gutzwiller projec- 
tion onto states with no double occupancies. As a result 
of the projection, 'F can be written as a sum over bond 
configuations of singlets and hole pairs which cover the 
lattice as depicted in Fig. 

u(jCij) and v(vij) are independent spin and hole bond 
parameters respectively. Uij > connects i on sublattice 
A to j £ B respectively, which ensures Marshall's sign. 
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FIG. 1. A bond configuration in the doped RVB states 
v]. Solid (empty) circles represent spins (holes) with 
bond correlations Uij (vki)- 
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The expectation value of an observable is computed 
by a sum over loop coverings: 

(O) = J2 w ^Oa^, 

7 ,A 7 



^A,7 = 



^-det 2 ^! TT 2 TT 



(2) 



AeA 7 \ (ij)GA 



A 7 denotes a list of directed loops {X a } which cover the 
lattice except for subset 7 of xL 2 hole sites. W are pos- 
itive Boltzmann weights, with which the Monte Carlo 
Metropolis step is defined following Refs. J^Q. Ergodic- 
ity and convergence of our pro gram was checked against 
precise transfer matrix results [ p.3| . 

The Gutzwillcr Approximation (GA) amounts to drop- 
ping the projector V(x) in state (||) and setting ^ ^ tp, 
after adjusting the overall normalization of u and v to 



satisfy the global constraints 



(1 — x), and 



(nj ) = x. ip is a Fock state of decoupled spins and holes 



V. 



with easily computable correlations |11 

For the spin correlations we use power law decaying 
functions u p (r) = l/r p . The single variational param- 
eter p determines the long range spin correlations and 
local magnetization (Other forms for u(r), with 

qualitatively similar results, will be described elsewhere 

We discuss four cases of inter-hole correlations: 



v ins( r i. 



Vmet(r) 

v a (r) 



G 7 
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(3) 



= s, d 



where fj are nearest neighbor vectors on the square lat- 
tice, c s = 1 and Cd = f] 2 — f) 2 . 

v lns P u f s t ne xL 2 holes on random sites. This state 
describes an insulator with disordered localized charges. 
i e t has filled Fermi pockets S, containing xL 2 occu- 



'V, 



pied k-states centered around (±7r/2, ±7r/2) in the Bril- 
louin zone pj|| . It describes weakly interacting holes in 
a "metallic" state. In real space, v me t{^ij) decays slowly 
as ~ r~ 3 / 2 . Correlations in this state were previously 
computed by Bonesteel and Wilkins 0. 

v s and Vd describe tightly bound hole pairs in relative 
s and d-wave symmetry respectively. 
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FIG. 2. (a) The local magnetization squared of doped RVB 
wavefunctions ^\u v , v] versus the variational power p, defined 
by the bond parameters Up(r) = l/r p . Lattice size is 40 x 40, 
and hole concentration is 10%. Results are agree well with the 
Gutzwiller approximation (solid line). The hole bond param- 
eters v are defined in Eqj^. Note that mo is weakly dependent 
on v; the data for v = v 3 , overlaps that of Vd- (b) Finite size 
scaling of toq(L) for p = 3.3 which indicates vanishing local 
magnetization at L — > 00. 

Order parameters: The local magnetization mo(L) is 
defined by m 2 = l/Z 2 £ i (S<S J -) I ,e~ i{ ' r,,r)ry - In Fi g|(a), 



ml(p) for 



0.1] is plotted for various choices of 



v. The GA (solid line) seems to work well for mo(p). Fi- 
nite size scaling in FigJ|(b) indicates vanishing long range 
order mo — > at p c — 3.3, which lowers the bound given 
previously by Ref. ||: at p c < 5. The GA approximation 
at L — > 00 suggests that p c < 3 |^2| . 

The superconducting singlet order parameters are 



Ay = flfkdibj - b iaj )/V2 



(4) 



By gauge invariance imposed by the Gutzwiller projector, 
(A Sjd ) = 0. However, 9[u, v s ( d )', x > 0] describe true s 
(d)-wave superconductors as seen by the (off-diagonal) 
long range order in A Sl <j flU^f . 

In contrast, the insulator states ^>[u,Vi ns ,x] and the 
"metallic" states ^>[u,v m et,x] have no long range super- 
conducting order of either symmetry Ji2fl . 

Effective Hamiltonians: Magnetic order is driven by 
the diluted Heisenberg model Q, 



u' ./^s..s,.i /;/, Ki /;/ ;! 

(ij) 



(5) 



where e.g. S" 21 + iS* = a^6. In Fig. || the expecta- 
tion value E mag (p) — (H J ) is plotted as a function of 
TOq(p) f° r x — 0-1 an d various choices of v from (^). 
Within numerical errors all states minimize 7i J at around 
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Vmin ~ 2.7, which by Fig. ||(a) yields local magnetiza- 
tion of tuq — 0.08 w mg(0)(l — 2x), where m (0) agrees 
with the ground state local magnetization of the undoped 
Heisenberg model We have found that p m in ~ 2.7 
appears to be independent of x in for < x < 0.15. 
Thus we conclude that aside from the trivial kincmatical 
constraints, the hole density and correlations have little 
effect on the magnetic energy at low doping. 

A single hole hopping in the antiferromagnetic back- 
ground has been shown by semiclassical arguments 
p7| , ^0| , to be effectively restricted at low energies to hop- 
ping between sites on the same sublattice: Next we con- 
sider the single hole hopping process 

W= E tkflfMai + blbi) (6) 

where i, k are removed by two adjacent lattice steps, and 
t' > 0. Unconstrained, the single hole ground state of 
if* has momentum on the edge of the magnetic Bril- 
louin zone, in agreement with exact diagonalization of 
t — J clusters Q. (Note: mier-sublattice hopping (the 
t-term in the t-J model) is a high energy process in the 
antiferromagnetic background. It primarily renormalizes 
t' and the quasiparticle's short distance structure |]l7| , |l0| . 

The single holes hopping (||) prefers the metallic state 
v = v met over the superconductor v — v s , v& (l2) . It 
also prefers longer range u(r) and thus actually enhances 
magnetic order at low doping. This is a type of a Nagaoka 
effect, where mobile holes seperately polarize each of the 
sublattices ferromagnetically. 

Now we consider Cooper pairs hopping terms 

^' = -J'K>j J .A i& + X a;a,] (7) 

The first term is derived from the large U Hubbard model 
to order J' = t 2 /U ll(J. It includes a rotation of the sin- 
glet pair, which prefers Vd over v s . The second term is 
a parallel translation of singlets. It prefers superconduc- 
tivity with v = Vd over metallic states with v = v me t 

& a n 

In Fig. |3| the ground state energy E p h of (|7]) is plotted 
for v — Vd, £=0.1 and L—A0. For v = v s , E p h > 0. 
The variational energy is minimized at p — 3.35, which 
by the finite size scaling of Fig.|2|(b) indicates vanishing 
too a t large L. Note the striking difference between the 
minima of E mag and E p h- Thus, Cooper pair hopping 
drives the groundstate toward a spin liquid phase! 
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FIG. 3. The magnetic energy E mag (y) and Cooper pair 
hopping energy E v u (^) versus local magnetization squared 
mo, used as a variational parameter. The density of holes 
is 0.1 and lattice size is L — 40. The magnetic energy is 
minimized at jiiq « 0.08, consistent with a diluted quan- 
tum Heisenberg antiferromagnet, and is weakly dependent 
on inter-hole correlations. The points of E mag (v s ) overlaps 
E mag (vd). In contrast, Cooper pair hopping prefers vanish- 
ing mo at L — > oo. 

A simple explanation is that pairs can hop with greater 
overlap when parallel bonds have maximum singlet com- 
ponents. When u(r) is longer ranged, triplet contribu- 
tions are larger, which inhibits pair derealization. In- 
cidentally the Gutzwiller approximation (GA) fails to 
predict this effect since it decouples the local correlations 
between spins and hole pairs. 

Since 7i J is the effective model which drives supercon- 
ductivity it produces phase stiffness, which in the contin- 
uum approximation is given by 

H J 'nY±J ^(V^) 2 (8) 

The stiffness constant Vq can be determined variation- 
ally from the RVB states, by imposing a uniform gauge 
field twist on the bond parameters Vij — > Vijexp(i(xi + 
Xj)4>/2L) and measuring E p h((j)) to find Vo = d 2 E p h / def) 2 . 

Following Ref . || , at low doping for the square lattice 
Vq is roughly equal to T c . 

results In Fig.(Q) we show our main result: The stag- 
gered magnetization too for H. J + H. J is plotted against 
the superconducting to magnetic stiffness ratio Vo/ J for 
different doping concentrations x = 0.05, 0.1, 0.15. 

Two primary observations are made: (i) The local 
magnetization is sharply reduced at relatively low su- 
perconducting stiffness (and T c /J). (ii) The relation be- 
tween too and Vq/ J appears to be weakly dependent on 
the precise hole concentration. 
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FIG. 4. The relation between thermodynamic local mag- 
netization rriQ =oc and superconducting phase stiffness Vb (re- 
lated to T c , see text). J is the Heisenberg exchange energy. 
The points are considered upper bounds on mo, which may 
even vanish for Vo/J > 0.2. 

Discussion: Because of finite size uncertainty, mo in 
Fig.(^) is an upper bound on the thermodynamic local 
magnetization. The GA extrapolation suggests that Too 
may actually vanish already Vq/J > 0.2. This is in qual- 
itative agreement with the doping dependent of the local 
magnetization measured by Refs. which diminishes 
rapidly above the onset of superconductivity. 

In a quantized theory of stripes [^0) , mechanisms for 
diminishing too assume anisotropic magnetic couplings, 
or fluctuating anti-phase domain walls. A direct connec- 
tion between superconductivity and mo is not obvious in 
these approaches. 

In a recent projected SO(5) theory spins and hole- 
pairs dynamics have been considered with excluded dou- 
ble occupancies. A variational relation is obtained be- 
tween superconducting stiffness and the magnetic order 
parameter, which resembles the results of this paper. 
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